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Abstract 



The excitation spectrum of the two-chain 5 = 1/2 Heisenberg spin lad- 
der with additional second neighbor frustrating interactions is studied by a 
variety of techniques. A description, based on a mapping of the model onto a 
Bose gas of hard-core triplets is used to determine the one- and two-particle 
excitation spectra. We find that low-lying singlet and triplet bound states 
are present and their binding energy increases with increasing frustration. In 
addition, many-particle bound states are found by exact diagonalization and 
variational methods. We prove that the larger the number of bound particles 
the larger the binding energy. Thus the excitation spectrum has a complex 
structure and consists of elementary triplets and composite many-particle sin- 
glet and triplet bound states. The composite excitations mix strongly with 
the elementary ones in the coupling regime where quantum fluctuations are 
strong. The quantum phase transition, known to take place in this model 
at critical frustration is interpreted as a condensation process of (infinitely) 
large many-particle bound states. 



Typeset using REVTpjX 



1 



I. INTRODUCTION 



The S — 1/2 quantum spin ladder is relevant to a number of quasi one-dimensional 
compounds^ and the list is growing as more materials become of experimental interest. 
Theoretically the two-leg ladder is, due to its geometry, the most simple realization of 
a "spin-liquid" - a quantum disordered state with gapped elementary excitations. The 
excitation spectrum of the ladder has been analyzed by a variety of techniques, including 
weak-coupling field theory mappings!, exact diagonalization of small clustersi'S and density- 
matrix renormalization group (DMRG) studies!. Also, strong-coupling techniques have been 
extensively used, such as dimer series expansions to high orders! and mapping onto effective 
bosonic theorieJBi'i. The dispersion of the lowest triplet excitation as well as the gap in the 
spectrum are quite well understood within the aforementioned approaches. 

Recently, an additional branch of excitations - two-magnon bound states were found in 
the spin ladder mo deli. Such bound states were also predicted for the dimerized quantum 
spin chain0 which is another quantum system with a disordered ground state and gapped 
excitations. Bound states in quasi one- dimensional gapped spin systems have also been 
observed experimentally^ although it is still not clear which one (or perhaps a combination 
of the two above@) is the relevant model for their description. Two of us have recently 
pointed ouM that bound states exist, in fact, in all one and two-dimensional quantum spin 
systems with dimerization of which the spin ladder and the dimerized chain are particular 
examples. 

In the present paper we study the two-leg spin ladder with additional second neighbor 
frustrating interactions between the chains. This model was introduced quite recently and 
analyzed numerically via dimer series expansions^, DMRG0 and exact diagonalizationslil'EiJ. 
A quantum phase transition was found as frustration increases from an antiferromagnetic 
(AF) ladder into Haldane (ferromagnetic ladder) phase. The excitation spectrum changes 
dramatically as one approaches the quantum transition pointlll. In a coupling region before 
the transition a singlet state appears in the triplet gap and at the transition both triplet and 
singlet gaps seem to approach zercll'0. We will show in the present work that as frustration 
increases a number of low-energy many-particle bound states appear in the spectrum which 
mix strongly with the one-particle excitations. The energies of the bound states decrease 
with increasing frustration and number of particles forming them. Thus the quantum tran- 
sition can be viewed as softening of a very complex excitation, composed of many-particle 
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bound states. 

Consider the Hamiltonian of two coupled 5 = 1/2 chains (spin ladder): 

# = E { -kSi-s; + j (s,.s m + s;.s^) + j 2 (s,.s: +1 + s:.s i+1 )} , (i) 

i 

where the intra-chain (J) and the inter-chain (Jj_, J 2 ) interactions are assumed antiferro- 
magnetic J, J 2 , J± > 0. In Eq.(l) J 2 is a second neighbor inter chain coupling which causes 
frustration. In order to analyze the excitation spectrum of (1) it is convenient to adopt the 
strong-coupling viewpoint. At large J± 3> J, J 2 the ground state consists of inter-chain spin 
singlets \GS) = |1,0)|2,0)|3,0)..., where \i, 0) = [| t)i| I)- - I !)•]• Since each singlet 
can be excited into a triplet state it is natural to introduce a creation operator V ai for this 
excitation: 

\h a ) = *LM)> a = x,y,z. (2) 

The representation of the spin operators in terms of t^ ai was introduced by Sachdev and 
BhattS: 

Si 2 = —(^a i tf a — ie a /3ryt ptry) . (3) 

After application of this transformation to (1), or, equivalently, after calculating the matrix 
elements of the "hopping" terms J and J 2 in (1), we find: 



H= £ \jj ai t 



2 

where we have defined 

A = J-J 2 , /i = J+J 2 . (5) 

In addition, we have to restrict the Hilbert space by introducing the following hard-core 
on-site constraintS 

tU% = 0. (6) 

This exclusion of double occupancy reflects the quantization of spin and ensures the unique- 
ness of the mapping from (1) to (4). 
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The Hamiltonian ([]]) as well as (|J) is symmetric under permutation of the ladder legs. 
Therefore all excitations can be classified according to this symmetry. Following standard 
notations we will denote the antisymmetric excitations (k± = ir) by the index u and sym- 
metric ones {k± = 0) by the index g. It is clear that the operator t a % (elementary triplet) 
corresponds to the w-excitation. 

The rest of the paper is organized as follows. In Section II we describe the one-particle 
(triplet) excitation spectrum. In Section III the two-particle problem is considered and 
bound states in various channels are analyzed. Section IV addresses the bound state problem 
for many particles focusing mainly on the case of three particles. Section V presents our 
analysis of the quantum phase transition in light of the previous results and summarizes the 
work. 



II. ELEMENTARY TRIPLET 

At the quadratic level the Hamiltonian (|j) can be diagonalized by a combination of 
Fourier and Bogoliubov transformations t a k = u^iak + ^fc^l_fc- This gives the excitation 
spectrum: oj\ = A\ — B\, where Ak = J± + X cos k and Bk = A cos k. We find, in agreement 
with previous work@i, that the effect of the quartic terms in (4) on the triplet spectrum 
is small and therefore we proceed by treating these terms in mean field theory. This is 
equivalent to taking into account only one-loop diagrams (first order in fi). These diagrams 
lead to the renormalization: 

Ak = J± + (A + 2/x/i) cos k, Bk = (A — 2figi) cos k, (7) 

where 

fx = (tUai+x) = N- 1 v\ cos q (8) 

q 

9x = (taitai+x) = A^^u^cosg. 

q 

The above corrections are numerically quite small. The dominant contribution to the 
spectrum renormalization is related to the hard core condition Eq.(6). This condition is 
typically taken into account in the mean-field approximationl'0. The latter is essentially 
uncontrolled, especially for a quasi- ID system. To deal with the constraint we will use 
the diagrammatic approach developed by us in Ref. [17]. An infinite on-site repulsion is 
introduced in this approach in order to forbid the double occupancy: 
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Hu = jY, t U%tf3itc tl , U^oo. (9) 

Since the interaction is infinite, the exact scattering amplitude r a p tJ g(K) = 
T(K) {S ai Sf3s + 5 a s5^), K = (k,oo), for the triplets has to be found. This quantity can 
be found by resumming the infinite series shown in Fig. 1(a). One can easily see that V 
depends on the total energy and momentum of the incoming particles K — K\ + Ki- The 
interaction (Q) is local and non-retarded which allows us to obtain the analytic expression0'i 
(in the limit U — > 00) 

d 2 Q 



(2*) 

= 1 U q U k-q | 1 V q V k-q /^n 

N q 00 — 00 q — 00 k - q N ^ 00 + 00 q + UJ k _ q ' 

Here G(Q) is the normal Green's function (GF) G(k,t) = -i(T(t ka (t)t\ a (0))): 

G{k,u,) = ^— (11) 

00 — 00 k + 10 00 + 00k — to 

and the Bogoliubov coefficients ul,vl = ±1/2 + A k j2oo k . The basic approximation made in 
the derivation of T(K) is the neglect of all anomalous scattering vertices, which are present 
in the theory due to the existence of anomalous GF's, GA(k,t) = —iiTit^^^tl^O))). 

n n \ u kVk u k v k 

G A (k,u) = — ■ (12) 

00 — oo k + 10 00 + oo k — w 

Our crucial observation!^ is that all anomalous contributions are suppressed by a small 
parameter which is present in the theory - the density of triplet excitations n t = Y^aVai^ai) — 
3iV _1 E g ^- We find that n t « 0.1 {J±/J = 2), n t w 0.25 {J ± /J = 1) and it generally 
increases as J± decreases. Since summation of ladders with anomalous GF's brings additional 
powers of v q into T, their contribution is small compared to the dominant one of Eq . (|1C|) . 
In the following analysis we will take into account only the contributions to the self-energy 
which are at most linear in the triplet density n t and therefore we also neglect the second 
term in Eq. flT0|) . Thus our approach is expected to work as long as the gas of triplets is 
dilute enough {n t is small). 

The normal self-energy which includes only the first power of the amplitude V is given 
by the diagram in Fig. 1(b): 

S (Br) (k, u) = 1 £ v 2 q T(k + q,oo-oo q ). (13) 

q 
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This is the dominant contribution to the spectrum renormalization as emphasized by 
BruecknerB who developed the technique described above in order to study systems of 
strongly interacting fermions. 

In the dilute gas approximation there are other diagrams which are formally at most 
linear in n t but still numerically give contributions much smaller than the one of Eq . (fL3|) . 
The first one is the "rainbow" correction to the anomalous self-energy which is proportional 
to y/nt and is shown in Fig. 2 (a): 

^ = 4EW(0,0). (14) 

q 

This anomalous self-energy enforces the condition 

(tltl) = N- 1 Y,u k v k = 0. (15) 

k 

The parameters Uk and v k found in the zeroth approximation do not satisfy this condition. 
Taking into account the self energy ([H]) gives the corrected values of Uk and Vk which do 
satisfy (|I3|). This can be seen from the formula for the renormalized Bogoliubov coefficients, 
Eq.(^3l) below. Since (|i~4|) is independent of k and u, technically one can take into account the 
anomalous self-energy by introducing the term A Yn, a (*L*L + tadcd) into the Hamiltonian 
(U) and choosing the Lagrange multiplier A from the condition (|i~5l). 

The next correction is the contribution to the normal self-energy given by the diagram 
shown in Fig.2(b), where the square denotes the scattering amplitude (|K]). A standard 
calculation gives the expression for this diagram 

N 2 p>q UJ-U p -UJ q - UJk+p-q 

Another correction is given by the diagram shown in Fig.2(c) plus the same diagram but 
with the positions of T and fi reversed. The result is 

E (2c) , k x = _^ ±_ y cos(p - q){u p v p ){u q v q )u 2 k+p _ q T(k + q,u- u q ) 

N 2 U-U p -UJ q - UJk+p-q 

The last correction linear in the triplet density is shown at Fig. 2(d). The corresponding 
expression is 

^ ( 2d)n \ 1 V- ( u p v p)( u q v <i) u k-p-i u k-q-i r ( k -P,u-UJp)r(k-q,u-UJ q ) / 1S rt 

N p~q,l {U — Up-Ul- Uk-p-l){uj -Ug — Ul- LOk-q-l) 

x [8 cos(p + q + I — k) + 10 cos(p — q)]. 
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Let us stress again that all normal self-energy contributions (|T3|) , (|T6|) , ([H]) , (|18D are quadratic 
in v g and hence linear in the triplet density. The anomalous self-energy (|T3|) is linear in v q 
and thus proportional to ^Jn~ t . 

In order to find the renormalized spectrum, one has to solve the set of two coupled Dyson 
equations for the normal and anomalous GF's, shown symbolically in Fig. 3. The result for 
the normal GF is: 

r(K] = u + A k + X(-K) 

1 ) [co + A k + E(-K)][u;-A k -E(K)} + [B k + E A (KW 1 ' 

After separating this equation into a quasiparticle contribution and incoherent background, 
wc find0: 

G(k, u) = Z f\ - Z f A + G - ( 2 °) 
00 — i l k + to UJ + \ l k — 10 

The renormalized triplet spectrum and the renormalization constant are: 



n k = Z k J[A k + 0)] 2 - [B k + £ A ] 2 , (21) 



u>=0 

Here the normal self-energy operator is given by Eqs . ( |I3[ ) , (|i~6] ) , (|17D , dl8|) 

E(k, uj) = £ (Br) + £ (2b) + S( 2c ) + (22) 



and the anomalous self-energy operator is given by Eq. (|14|) . The renormalized Bogoliubov 
coefficients in (^Op are: 



ulv , = A + ^t±pm . ( 23 ) 

k k 2 2Q fe 



Equations fll0|,|2"2l , [2l| , |23D have to be solved self-consistently for E(/c, 0) and Z k . From Eq.(pO| 



it also follows that one has to replace u k — > y/ZkUk, v k — > ^fZ^Yk in all expressions presented 
above and below. 

Let us demonstrate how this approach works in the strong-coupling limit Jj_ ^> /i, A. To 
first order in \/J±, A k = J± + Acosfc and B k = Xcosk. This leads to u k ~ A k , u k ~ 1, 
v k « — (A/2Jj_) cosfc and fi — 0, gi — — X/AJ±. Substitution into (PH|), ([HJ), and ( |2"2"| ) gives 

r(fc,o;)= 2J ± -w, (24) 
E(ife,a;)= £ (Br) (A;,cu) = ^(\/ J±) 2 (3J ± - oj). 



Note that self-energy corrections (0), (|16"D, (|T7|), and (|1|) do not contribute in this order. 
Then from Eq.(|2~I|) we find the quasiparticle residue Z = 1 — (1/2)(A/ Jj_) 2 and the dispersion 

3A 2 A 2 

&k = J± + A cos k + — — — cos 2k. (25) 

4jj_ 4J_|_ 

The result (|25 ) agrees with that obtained by direct 1/ Jj_ expansion!!! to this order. 

It is also useful to consider the next order in 1/J±. Using the first order calculation 
presented above we find A k = J± + A cos k and B k = A(l + fi/2J±) cos k and hence u k ~ 1, 
fit ~ —B k /2A k « — (A/2Jj_)(l + /i/2 J_l)(1 — A/ Jj_ cos fc) cos A;. The scattering amplitude T is 
not changed in this order and thus given by Eq.(|24|). The anomalous self energy calculated 
according to Eq . (|14"D and the contributions to the normal self-energy given by Eqs.([ilf),(|l6l), 
(0), © are: 

E A = A 2 /2J ± (l + /i/2J ± ), 

S^) = i(A/J ± ) 2 (l + /U /J ± )(3J ± -c), 

S (2b) = ^cosA;, (26) 

yj(2c) Ai-^ 2 



E (2d) 



4J 2 ' 
5/iA 2 



8J : 



2 



Substituting these into Eqs.(^),(^l|) we find the elementary triplet dispersion to order 1/ J\: 

i x2 f 3 1 n1 \ A 3 / 1 ,1 A A 2 /i /3 1 

i Zfc = Jj_ + A cos k H cos 2k H ^ — cos k H — cos 3k H =- cos 2k 

J± V4 4 / Jl V 4 8 / JJ V8 4 

(27) 



Using Eq.(^) one can also prove that the condition (|15D is satisfied. The result (P7|) agrees 
with that obtained by direct 1/J± expansioni'lli to this order. 

The technique presented above is certainly not the simplest way to construct the 1/ J± 
expansion. Moreover it can not reproduce terms of order 1/Jj_ and higher because contri- 
butions to the self energies which are quadratic and higher order in the triplet density have 
been neglected in our approach. However the advantage of the method comes from the fact 
that n t remains relatively small (0.25) even for J/J± = 1. The purpose of the presented ex- 
ercise was to demonstrate that the result of our approach coincides with the result obtained 
by perturbation theory around the dimer limit to the relevant order. 
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For arbitrary J± a self-consistent numerical solution of Eqs. ( |10| , 13| , ^T , |23 ) is required. 
The triplet excitation spectra obtained from this solution for J±/J = 2 are shown in Fig.4. 
For comparison we present the spectrum for J 2 = which only includes the Brueckner 
correction (|13|) as well as the spectrum which includes all terms linear in nt (the self-energies 
(|I^),(|l^),(|l7D,(|l8D, in addition to (|T3l)). One can see that the Brueckner diagram is the 
most important one. All other corrections are much less important, however we will keep 
them in all subsequent calculations. Notice that the correlation corrections described above 
renormalize the spectrum very strongly as can be seen by comparing with the bare dispersion 
(all correlations neglected, U — /j, — 0): ojf. = J\ + 2XJ±cosk. The bare spectrum even 
becomes unstable for J± < 2A. In Fig.4 we also present for comparison dispersions obtained 
by 8-th order dimer series expansion!. The agreement between our calculation and these 
curves is excellent which reflects the smallness of the triplet density n t ~ 0.1. In Fig. 5 we 
present similar plots for the case J± = J. Looking at the curves at J2 = one can say 
that the agreement between our theory and the result obtained by series expansions is still 
reasonable because the triplet density in this case is n t w 0.25 and hence one has to expect 
about 25% disagreement. However as J2 increases the disagreement increases (especially at 
the point k = 0) in spite of the fact that according to our calculation the triplet density 
does not increase and even slightly decreases. Moreover the excitation energy at k = 
vanishes at J2 ~ 0.6J, which signals a quantum phase transition into the Haldane phase. 
Our calculation however does not give any indication of the triplet mode becoming soft at 
k = 0. Therefore something important is missing in our approach. We will demonstrate 
in Section IV that what is missing is the contribution of low-energy many-particle bound 
states (3,5,7... particles) which have w-symmetry and therefore can mix with the elementary 
triplet. 

Next, we proceed with the analysis of two-particle bound states which have ^-symmetry 
and therefore do not mix with the elementary triplet. 



III. TWO-PARTICLE BOUND STATES 

The quartic interaction in the Hamiltonian (^) leads to attraction between two triplet 
excitations. We will show that the attraction is strong enough to form a singlet (S=0) and 
a triplet (S=l) bound state. The method we employ essentially follows our previous worki. 

Consider the scattering of two triplets: q±a + q2f3 — > q?^ + q^S and introduce the total 
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iff) and relative (q) momentum of the pair q\ = Q/2 + q, q 2 = Q/2 — q, g 3 = Q/2 + p, and 
Q4 — Q/2 ~ P- The bare (Born) scattering amplitude is (see Fig. 6(a)): 



M at3n s = yU (Sajfyd - Sa^S-ys) COs(q + p) + 

/j, (SagSfr - 5 a p5 lS ) cos(g - p) + 

U (SajSpg + 5 aS 5pj). (28) 

The fi and the U terms arise from the quartic interaction in (Q) and the constraint @ 
respectively. We also have to take into account that the triplet excitation differs from the 
bare one due to the Bogoliubov transformation and the quasiparticle residue. Therefore the 
following substitution has to be made: 

M a/3 ^s — » \l^iU qi ^^ 2 U q2 ^J^JJ q3 ^J^JJ qA M a p^ s . (29) 

The bound state satisfies the Bethe-Salpeter equation for the poles of the exact scattering 
amplitude M. This equation is presented graphically in Fig.6(b) and has the formll: 



En — fio/oin — Vl 



Q/2+q ~ "Q/2-g 



^(q) = lJ^M(Q,q jP )^(p). (30) 



Here M(Q, q,p) is the scattering amplitude in the appropriate channel, Eq is the energy of 
the bound state and ip(q) is the two-particle wave function. The factor of 2 in Eq.(|3TJ|) is 
related to the symmetry of the diagram on the right hand side of Fig. 6(b) under the exchange 
of the two intermediate lines. Thus in order to avoid double counting of the intermediate 
states, the result has to be divided by two. Let us introduce the minimum energy for 
two excitations with given total momentum (lower edge of the two-particle continuum) 
Eq = mhig + Qg/2-gj. If a bound state exists then its energy is lower than the 

continuum Eq < Eq. The binding energy is defined as eg = Eq — Eq > 0. 
In the singlet (S=0) channel the scattering amplitude is: 

M (0) = -5 af3 5^M al3r/S = -Afi cos qcosp + 2U. (31) 
3 

First, consider the strong-coupling limit J± ^> J, J 2- Let us keep terms up to first order in 
1/ J±, i.e. take Q q from Eq.(f23|). The lower edge of the continuum in this order is: 



E h = 2J 1 + ^ + -^cosQ-2Acosg/2 ,Q<Q* ^ 

1J: - 1 +2^cosQ + J ± (cos 2 Q/2)/cosQ ,Q > Q* 
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Here Q* is determined from the equation: (cosQ* /2)/ cosQ* = —\/J±. Notice that in the 
strict limit X/J± = one has Q* = tt and thus the upper line in Eg. (P^|) is sufficient. The 
equation for the bound state reads: 



3A 2 A 2 

Eq^ — 2 J_l — 2 A cos Q/2 cos q 1 cos Q cos 2q 

2J i 2J 



ip(q) 

dp [dp 



f dp [dp 
= — 2/icosg / — cos pip(p) + U / — ip{p)- (33) 
J 2n J 2ir 

Since we work to order l/J± and both Z q , U q — 1 + 0(A 2 /J 2 ), these quantities have been 

set to unity in fl3H| ). Due to the infinite repulsion (U — > oo), a Lagrange multiplier has to be 

introduced to enforce the condition / dpip(p) = (meaning that the bound state is d-wave 

like). The solution of Eq. fl3"3]) to leading order for the wave- function and next to leading 

order for the energy is: 

^( ftq )-V^ IT »l±gL_ + (^) ,34, 

Ef = 2J ± + ^- - Ml + C 2 Q ) - ^-(1 + CD cosQ (35) 

where we have introduced the notation 

C Q = -cosQ/2. (36) 
/i 

Thus we see that in the strong-coupling limit a singlet bound state always exists. At J± = 2 J, 
J2 = Eq.flHOp with the substitution ( |2"§D has to be solved numerically and the result is 
presented in Fig. 7. We find that for k < 2tt/5 the binding energy is practically zero in this 
case. 

In the triplet (S=l) channel the scattering amplitude is: 

M (1) = ^e pa pe pl sM al3 ^s = -2fi sin q sin p. (37) 

In this formula there is no summation over the index p which gives the spin of the bound 
state. By solving Eq.(^) in the limit Jj_ ^> J, J2 we obtain for the wave-function and the 
binding energy: 

V><% Q) = v/l/2-2C 2 f 9 + O (^-) (38) 

v v 1/2 + 2Cq + 2Cq cos q \ fj,J± J 



^ 1) =2J ± + ^-^(l+4C 2 )-^(6C 2 -l/2)cosQ ! C Q < 1/2. (39) 
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For Cq > 1/2 we find that the binding energy vanishes, e^q = Eq — Eq^ = 0, which means 
that at J 2 = (/i = A = J) the triplet bound state only exists for momenta k > Q c = 2tx/2> 
(in the strong-coupling limit)!. At J± = 2 J, J 2 = the numerical solution of Eq. (|30|) , plotted 



in Fig. 7 (with the additional contribution Eq.(|4T|)) shows that the bound state exists down 
to k w 7r/2. 

Finally, we find that there is no bound state in the tensor (S=2) channel. This is due to 
the fact that the scattering amplitude in this case = 2ficosqcosp + 2U corresponds to 
repulsion and consequently there is no solution of the Bethe-Salpeter equation with positive 
binding energy. However a solution exists with energy above the upper edge of the two- 
particle continuum. In the simplest case J = J2, X = we find to leading order E^ = 
2J_|_ + /i/2 and thus the "anti-binding" energy is /x/2. 

Equation (|3(]) takes into account the potential interaction between two dressed elemen- 
tary triplets, but it does not take into account the contribution of quantum fluctuations into 
binding. Let us consider this effect. In the strong coupling limit the first correction to the 
ground state energy of the system is due to the term li^L+i m ^ ne Hamiltonian (f|) which 
virtually excites a pair of triplets. Thus the energy correction per link to lowest order is 

SEo = -3^2, (40) 

where the coefficient 3 is due to the number of possible polarizations^. When we have 
a state with a real elementary triplet, the quasiparticle (triplet) blocks virtual excitations 
on two links and this increases its energy by 2|5J5o|- This is the physical origin of the 
third term in the dispersion (p5f). Now let us consider two quasiparticles. When they are 
separated by more than one lattice spacing they block four links, but when they are on 
nearest neighbor sites they block only three links. This gives an effective attraction 5Eq. 
However two quasiparticles in a singlet (S=0) state can virtually annihilate because of the 
term in the Hamiltonian (|j) which has the same tensor structure. This term gives 

—SE and consequently the net effective attraction due to quantum fluctuations vanishes. 
For the triplet (S=l) bound state there is no annihilation and therefore the energy level 
shift due to blocking of quantum fluctuations is@ 



dE 1 }? = SEn 



I vising (?, Q)^ 



2 



(41) 



The integral gives the probability amplitude for two quasiparticles to be on nearest neighbor 
sites. The two-particle triplet (S=l) bound state energy for J± = 2J, J 2 = is plotted in 
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Fig. 7 where the potential contribution as well as Eq. (f4T|) have been taken into account. 
While in the strong coupling limit the binding in the triplet channel is weaker than the one 
in the singlet channel (as can be seen from Eqs.fl3"5|),(|3"PD), for J± = 2 J, J 2 = the additional 
attraction due to blocking of quantum fluctuations pushes the triplet below the singlet for 
the range of momenta 47r/5 < q < ir. 

The sizes of the bound states can be determined from the corresponding wave functions. 
As expected the size increases with decreasing binding energy and near the threshold we find 
Rrms ~ ( e ) 2 > e ~~ * 0. The self-consistent evaluation of the sizes shows that both bound 
states typically extend over a few lattice spacingsi. 

The quantity which is directly measurable in inelastic neutron scattering experiments is 
the dynamical structure factor: 

S g , u (k,u) = J e^(Sr(k,t)SI' u (-k,0))dt, S% = S z>i ± S' Zji (42) 

The superscript corresponds to transverse (along the rungs) momentum k± = 0, n, i.e. 
S z 'i = S Zi i ± S' zi . The symmetric combination (k± = 0) gives the magnetic moment of 
the elementary triplet which is equal to unity. Therefore expressed in terms of Cartesian 
components the magnetic moment has the form M M = —ie^pv^p. This immediately gives 

S z ,i + S' z ,i = ~^ e zaptl li ti3i — > —ie za /3 ^ ^q^fc-VL^k-fp (43) 

1 

where we also have taken into account the Bogoliubov transformation. By projecting this 
operator onto the bound state wave function we find the contribution of the S = 1 bound 
state to the static structure factor S g (k) = J S g (k,uj)dLj/2TT: 

n 2 



S g (k) 



— ^ (q,k)u k /2+qVk/2- 
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~(A/ J ± ) 2 sin 2 k/2(l - ACl) + O (A 4 /Jl) • (44) 



In this formula Ck is defined by Eq . (|36|) . The substitution (uk,Vk) —> y/Zk(Uk,Vk) has to 
be made according to ( p0|) in order to find the result for arbitrary J±/J. We have also 
presented the leading order of the strong coupling expansion. 

A similar calculation in the u channel, i.e. for the elementary triplet gives 

S u {k) = K + v k f = 1 - A cos k + O (A 2 / jQ . (45) 

For J 2 = 0, J± = 2 J we have found by numerical evaluation of the corresponding expressions 
that S g (ir) / S u (tt) ~ 0.05 and thus the experimental signal is expected to be about 20 times 
weaker for the bound state compared to the elementary triplet!. 
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IV. MANY-PARTICLE BOUND STATES 



Let us first consider a three-particle bound state with total spin S=l (triplet). This state 
consists of an odd number of elementary triplets and hence has w-symmetry. A convenient 
way to solve the three-particle problem is to use the variational method. First consider the 
simplest ansatz: three triplet excitations on nearest neighbor sites. Such ansatz is valid in 
the limit of zero hopping (A = 0). A straightforward minimization of the expectation value 
of the Hamiltonian @ gives the energy and the wave function of this state: 



p(k\H\k) p = 3J ± - 1.25/x, 

where k and p are the momentum and the polarization of the state. Next, one can extend 
this ansatz by allowing each triplet to hop onto a nearby site (first order in A): 

4> p (k) = a\k) p + b\k)' p , (47) 
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The state ip p {k) must also be normalized, i.e. a 2 + b 2 = 1. The Hamiltonian has to be 
calculated in this basis, and additionally the energy level shifts due to blocking of quantum 
fluctuations have to be included, similarly to the discussion in the previous section. The 
result for the effective Hamilton matrix is: 



(H) 



eff 



^3^-1.25^ + 2^ £ ^ 

v 72 3J ± -p + lcosk+ 1 -l^ J 



(48) 



Notice that the quantum fluctuation correction in the second diagonal term ^ s slightly 

larger than the one in the first term. This is the same effect as the one discussed in the 
previous section - effective attraction due to suppression of quantum fluctuations. In this 
situation numerically this attraction is not very important. The energy of the three-particle 
bound state is 



hj-i(k) = SJt uH — cos AH \ i — i — wsm i -: 

SK ' 8 P 4 16 J ± \ \S 4 16 J ± J 2 



- + -cosA;H — H . (49) 



Consider first the strong coupling limit, J± ^> J, J 2 . For J 2 = (i.e. fi = A = J) eq. 
gives 
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E 3 (k = 0) = 3J ± - 1.68J, (50) 
E 3 (k = tt) = 3J ± - 2.09J. 

The state with = n is unstable with respect to decay into three elementary triplets because 
the energy of the elementary triplet is Q q = J± + Jcosq. However the state with k = 
is stable with respect to this decay. Nevertheless this state is also unstable since it can 
decay into a two-triplet bound state (Section III) and an elementary triplet. The threshold 
for this decay is 3J± — 2 J which is pretty close to E 3 (0) given by (|50|). Therefore a quite 
natural question arises: can improvements of the variational wave function push the energy 
E 3 (0) below the threshold? To check this we extended the ansatz (pE7|) by including states 
with double hopping (order A 2 ): t^^t^^, tJj_ 1 t^t^ +3 , and in-2*nVt-2- We find that E 3 (0) 
decreases to the value 3Jj_ — 1.77 J, but still remains above the decay threshold. Therefore 
we believe that in the strong coupling limit for J 2 = the three-particle bound state does 
not exist. However when J 2 > (0.3 — 0.4) J the bound state at k = becomes stable which 
follows immediately from Eq. ([49|) . 

For intermediate values of J± the three-particle state becomes stable for any J 2 . Let us 
consider three cases for J_l =2.1 . According to Eq.(^) 



(51) 



In all these cases any decay of the k = state is kinematically forbidden (this can be found 
from comparison with the elementary triplet and two-particle bound state spectra presented 
in Figs. 4, 7). 

Next, we compare the variational results with numerical exact diagonalization 
results we have obtained for a 2 x 10 ladder. Plots of the spectral function 
A(k, u) = — 7r _1 ImG'(fc, u + i5) in the w-channel (odd number of particles) for k = found 
by Lanczos diagonalization of the Hamiltonian (|1|) are presented in Fig.8. The first peak 
corresponds to the elementary triplet and the second one to the three-particle bound state. 
The positions of the second peak agree very well with Eq.(^TJ). For k = it we find numeri- 
cally that a second peak is absent for J 2 = 0, 0.4J whereas a peak with an extremely small 
spectral weight seems to exist for J 2 = 0.8J. This can be understood from the variational 
treatment since the state k = n can decay into three elementary triplets (compare (|51|) and 



E 3 (k = 


0) 


= 5.3J, 


E 3 (k = n) 


= 4.9 J, 


E 3 (k = 


0) 


= 4.4J, 


E 3 (k = n) 


= 4.2 J, 


E 3 (k = 


0) 


= 3.7J, 


E 3 (k = n) 


= 3.7J. 
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Fig. 4) for J 2 = 0. Even though this state is slightly below the threshold for J 2 = 0.4J, due 
to the limited accuracy of our calculation it is really hard to say whether it decays or not. 
However for J 2 = 0.8J the state k = n is well below the decay threshold, and indeed a peak 
exists in the corresponding spectral function. Thus we believe that the variational method 
captures quite accurately the main features of the spectrum. 

For Jj_ = J according to Eq.(fj9|) the three-particle bound state energy is 



J 2 = 0: E 3 (k = 


0) 


= 3.3J, 


E 3 (k = n) 


= 3.0J, 


J 2 = 0.4 J : E 3 (k = 


0) 


= 1.8 J, 


E 3 (k = n) 


= 1.6J, 


J 2 = 0.6 J : E 3 (k = 


0) 


= 1.2J, 


E 3 (k = n) 


= 1.1 J, 



Comparing with the exact diagonalization spectra presented in Fig. 9 and Fig. 10 one can see 
that the overall agreement is good. Notice that while for J 2 = the variational energies are 
higher than the numerical ones (as one would expect), for J 2 = 0.4 J, 0.6 J they are in fact 
lower. We attribute this effect to the mixing between the three-particle and the elementary 
triplet which has not been taken into account in our approach (see the discussion below). 

In the numerical spectra in Fig. 9 and Fig. 10 a third peak is also clearly seen. This 
is the five-particle bound state. To estimate its energy as well as the energies of bound 
states containing higher number of particles we could use the N = 00 approximation (N is 
the number of particles). In the limit A = the quartic term in the Hamiltonian Eq.([|) is 
identical to the Hamiltonian of an S — 1 Heisenberg chain with antiferromagnetic interaction 
fi/2. The ground state energy of the latter (for an infinite chain) is known quite accurately 
to be —0.700742^ per linkH. Therefore a crude estimate for the energy of an iV-particle 
bound state (containing N — 1 links) is 

E N = NJ ± - (N - 1) x 0.7//, (53) 

For the five-particle bound state by using the above formula and taking also into account the 
increase in energy due to blocking of quantum fluctuations (3\ 2 /J±), we obtain E 5 w 4. 5 J 
for J± = J, J 2 = 0, and E 5 ~ 1.9 J for J± = J, J 2 = 0.4 J, in qualitative agreement with the 
numerical results presented in Figs. 9, 10. 

Now we can address the problem formulated at the end of Sec. II: Why the diagrammatic 
approach developed in Sec. II, which works quite well for J 2 = 0, does not describe even 
qualitatively the triplet energy spectrum for J± = J and J 2 > 0?. In light of the results 
of the present section, we find that the essence of the problem is in the neglect of bound 
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states of three, five, etc. quasiparticles whose energies decrease with increasing J 2 . Indeed, 
let us fix J 2 = 0.4J and compare the energy of the elementary triplet at zero momentum 
from Fig. 5 (dashed line), f2 « 1.73J, with the energies of the three- and five-particle bound 
states E 3 (k = 0) ~ 1.8J, E 5 (k = 0) ~ 1.9J. They are quite close, and since all these 
states have the same quantum numbers they mix strongly. Notice that in the calculation 
of the one-particle properties as well as the three-particle problem we have not taken the 
mixing into account. Thus we expect the wave function in the w-sector (and similarly for 
the ^-sector) to be a superposition of states with different numbers of quasiparticles: 

= + Z 3 |^) (3) + Z 5 |^) {5) + . . . (54) 

In this situation the classification of the states by the number of "elementary" quasi particles 
is becoming meaningless, and the average number of excited triplets in the lowest excitation 
at k = is increasing. The full description of the energy spectrum requires the determination 
of the mixing coefficients in Eq.fl54|) which is beyond the scope of the present work and will 
be reported in the future. We expect that the energy of the "elementary" triplet will lower 
substantially at k = (with respect to the "naive" calculation of Sec. II) due to repulsion 
from the nearby many-particle bound states. In addition, as can be seen from the analysis 
of the three- and five-particle bound states, the larger J 2 the larger the number of many- 
particle bound states which have low energies and mix with the "elementary" triplet. In 
fact it becomes energetically more and more favorable to form states with larger and larger 
number of quasiparticles in them as J 2 increases. Thus we expect that the quasiparticle 
residue will decrease with increasing J 2 - an effect which indeed can be seen from our 
numerical analysis (see Fig. 9 for J 2 = 0.6 J). Eventually a situation may occur when the 
quasiparticle residue has vanished completely which means that very large size bound states 
completely dominate in the wave function Eq . (|54]) . This is the point where there is an 
excited triplet on every site and the ground state changes its nature. 



V. QUANTUM PHASE TRANSITION IN THE MODEL. SUMMARY AND 

CONCLUSIONS. 

The analysis of the previous section allows us to shed new light onto the nature of the 
quantum phase transition which takes place in the frustrated ladder model. The phase 
diagram of the model was determined in Ref.[13] and is presented in Fig. 11. At a critical 
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coupling J2c the ground state changes from that of an antiferro magnetic (AF) spin 
ladder to a ladder with an effective ferromagnetic interaction on the rungs (Haldane phase). 
From the point of view of the triplet excitations in the AF ladder phase, the Haldane phase 
is characterized by an excited triplet on every rung. Thus it is not surprising that bound 
states of many-particles become favorable energetically near the quantum transition point. 

The analysis of the energy spectrum is particularly simple on the line J2 = J(X = 0) 
where quantum fluctuations are absent completely. It is known that on this line there is 
an exact eigenstate of the Hamiltonian (H) which is a product of singlets (dimers) on each 
rung@. This is obvious from Eq.(H). This state is the ground state in the region J± > 1.4 J 
(see below). As J± decreases from a large value and approaches the quantum critical point, 
a number of singlet states appear in the triplet gap. Figure 12 presents a plot of the 
elementary triplet (ul), two-particle singlet (g2), three-particle triplet (u3) and four-particle 
singlet (g4). The energies of these states have been found by analytical diagonalization of 
the Hamiltonian Eq.(|): E ul = J ± , E g2 = 2J ± - 2 J, E u3 = 3J ± - 2.5J, E gA = 4J ± - 4.46J. 
It is clear that at the point J 2 = 2J the two-particle singlet crosses the one-particle triplet 
and thus becomes the lowest excitation in the system. Also we observe that the larger the 
number of bound particles the larger the rate of decrease of their energy. For comparison 
we have also schematically plotted the states u9 and glO. Thus we see that a number of 
singlets appear in the triplet gap and many level crossings take place. Notice that there is 
no mixing between the states since quantum fluctuations are absent (A = 0). At the point 
J± c = 1.4J the energy of the singlet composed of infinitely many quasiparticles becomes 
zero, E goo = 0, as can be seen from Eq.fl53|). The triplet (u) bound state energies do not 
cross the elementary triplet for any finite number of particles in them, however the infinite 
particle triplet becomes degenerate with the corresponding singlet E uoo = E goo = at the 
transition point (this also follows from Eq.(|53|)). 

We believe that the picture of the quantum transition presented above remains valid 
along the whole critical line (Fig. 11). The transition is characterized by softening of the 
singlet and triplet (at k = 0) modes which are basically very large size bound states of 
many quasiparticles in the appropriate channel. Slightly away from the critical line (on 
the AF side) the excitation wave function is a mixture of bound states with different num- 
ber of particles and the weight of the large-size bound states increases as the transition is 
approached. 

In summary, we have analyzed the properties of many-particle bound states in the frus- 
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trated ladder model. We have found that the excitation spectrum is quite complex and 
many-particle bound states are always present in the model. Frustration pushes the bound 
states to lower energies and the effective triplet and singlet spectra are very strongly renor- 
malized with respect to the simple ladder (no frustration). Thus the model is an ideal 
playground for studying complex excitations in quantum spin systems. 
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Equation ( |30|) can also be derived as a two-particle Schrodinger equation in momentum 
space H^ Q = E Q ^ Q . For the singlet bound state one has ^>q = J2 g ^°\q, Q)ta, qi ta, q2 \fy , 
where qi^ are defined in the text before Eq. (|28|) . For a triplet bound state with polariza- 
tion a: ^Q, a = Ytqip^Hq, Q) e a/3'yt Jf /3 )qi t\ i q 2 \0) ■ While this way of deriving the Bethe-Salpeter 
equation is quite simple and physically transparent, the approach presented in the text is 
more general and we give it preference since it is better suited for treating more compli- 
cated interaction terms (e.g. three-body interactions which appear in certain models with 
dimerization0). 

Strictly speaking Eq. ( [4"0"D is valid for J_l 3> J, J 2 . However it works quite well even for 
J± ~ J. For example at the point J± = J, J2 = 0, according to Eq.(|40|) 5E Q = —0.375 J 
while the numerical resulti'0 is 5E Q = —0.406 J. 

This effect is similar to the shift of energy levels due to radiative corrections in relativistic 
atomic physics (Lamb shift). 
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FIGURES 

FIG. 1. (a) Resummation of the infinite ladder for the scattering amplitude T. The dashed line 
represents the (infinite) two-particle interaction U. (b) The self-energy, corresponding to T. 

FIG. 2. Diagrams for the self-energy which contribute to linear order in the triplet density 
rif. The boxes represent the scattering amplitude T from Fig. 1(a). The wavy line stands for 
the two-particle interaction /u, Eq.(^|). Lines with a single arrow represent normal Green's func- 
tions (Eq.([ll|)) while lines with oppositely pointing arrows represent anomalous Green's functions 

(Eq.©)- 

FIG. 3. The coupled set of Dyson's equations for the normal and anomalous Green's functions. 
The anomalous self-energy (Fig. 2(a)) is denoted by A. The thin lines represent the bare Green's 
functions, Eq.([Tl]) (single arrow) and Eq.fll^) (double arrows). 

FIG. 4. The one-particle (triplet) excitation spectrum of the ladder for J± = 2 J. The solid dots 
represent numerical results obtained by dimer series expansions! for J2 = 0. The solid and dashed 
line are the results of the self-consistent numerical evaluation of the spectrum Eq.([2l|) for J2 = 
and 0.4J, respectively. The dotted line is the J2 = result when only the Brueckner self-energy 
Eq.(li~3l) is taken into account. 

FIG. 5. One-particle spectra for J± = J. The solid dots, open circles and solid squares are 
the dimer series expansion results of Ref.[13] for J2 = 0, 0.4 J and 0.6 J, respectively. The solid and 
dashed line are the results of the self-consistent numerical evaluation of the spectrum Eq.(|2l|) for 
J2 = and 0.4J, respectively. 

FIG. 6. (a) The bare (Born) scattering amplitude M. (b) the Bethe-Salpeter equation for the 
poles of the exact scattering amplitude M. 

FIG. 7. The excitation spectrum for J± = 2 J, J2 = including the singlet bound state (long 
dashed line) and the triplet bound state (dot-dashed line). The solid line E% is the lower edge of 
the two-particle continuum. 
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FIG. 8. Spectral function A(k, lo) for k = 0, J± = 2 J and several values of J2 obtained by 
Lanczos diagonalization of a 2 x 10 ladder. <5-functions are replaced by Lorentzians of width 0.1J. 

FIG. 9. Same as Fig.8 for k = 0, J ± = J. 

FIG. 10. Same as Fig.9 for k = tt, J ± = J. 

FIG. 11. Phase diagram of the frustrated ladder from Ref.[13]. The crosses represent the line 
J 2 = J where the ground state is a product of rung singlets. 

FIG. 12. Schematic excitation spectrum on the line J2 = J. 



23 



(a) 




FIG.l. 



+ 




FIG.2. 




FIG.3. 




FIG. 4. 




FIG. 5. 




FIG. 7. 



z — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — z 
- k - 




FIG. 8. 




FIG. 9. 




FIG. 10. 




FIG. 11. 




FIG. 12. 



